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Abstract 

In this letter we discuss the masses and the splittings of states in the spectrum 

of D and B mesons, as they appear in the approach of the QCD string with quarks at 
the ends. We find good agreement of our predictions with those of other QCD-motivated 
models as well as with the lattice and experimental data, including recent experimental 
results. We discuss the ordering pattern for P levels in D- and i?-mesonic spectrum. 

Data on the spectroscopy of heavy-light mesons coming from various experimental collab¬ 
orations are challenge for theorists, and these are D and B mesons to play an important role 
in checks of the validity and accuracy of the models. 

In this letter we address questions concerning the masses of orbitally excited D and B 
mesons in the method of the QCD string with quarks at the ends, paying special attention 
to the ordering of the P levels. The choice of D and B mesons is not accidental and is 
stipulated by the recent data on the masses and decays of the above mentioned heavy-light 
mesons coming from various experimental collaborations. Despite of the fact that there is no 
agreement between them and some resonances are not yet conhrmed, still we hnd it interesting 
to compare these experimental data, as well as those provided by other models and lattice 
simulations, with the predictions of our approach. First, let us remind the reader the basic 
ideas of the latter. 

Starting from the gauge-invariant wave function of the qq meson, 

^qq{x, y\A) = 'Fq-(a:)$(a:, y)'^q{y), (1) 

with $ being the parallel transporter, we write the Green’s function of the meson, 

^qq~ ^ 2 /( 2 ) 
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Meson 

a. Gel/" 

(Xg 

Co, MeV 

iriQ, MeV 

nig, MeV 

D 

0.17 

0.4 

196 

1400 

9 

B 

0.17 

0.39 

169 

4800 
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Table 1: Parameters of the Hamiltonian (|^)-(|^). 


and perform the integration over the quark and the gluonic helds. For the latter case we make 
use of the minimal area law asymptotic for the Wilson loop bounded by the quark and the 
antiquark trajectories (see, e.g., |Q), 


(ttP exp (^ig j>^ dz^A^ ~ exp (-aSmin), 


(3) 


where a is the string tension in the fundamental representation of the SU(3) colour group, and 
the area S'min can be approximated by means of the straight-line anzatz |^, 


^ ■ — 
*^min 


dt [ dpJ{liiw')^ — 

Jo * 

j3) = I3xi^{t) + (1 - I3)x2^,{t), 


(4) 


with Xi ^2 being the coordinates of the quark and the antiquark. Now, applying the Feynman- 
Schwinger representation to the single-quark propagators and introducing the einbein helds 
/ii ^2 to simplify the relativistic kinematics 0, we, hnally, arrive at the following expression for 
the Hamiltonian of the meson 


H — Hq + Vstr + Vsd, 


(5) 
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where in we supply the purely nonperturbative interaction, coming from the string-like 
picture of conhnement, by the perturbative Coulomb interaction (k = |q!s), as well as by 
the constant negative shift, Co, due to the light-quark self-energy P] strongly needed to bring 
the Regge trajectory intercepts into their experimental values. The term Vgtr deserves special 
attention, since it is originated from the square root in (§) and describes the contribution of the 
QCD string into the total inertia of the rotating qq system. This contribution is important to 
establish the correct slope of the mesonic Regge trajectories [0. We keep the hrst two terms in 
its expansion in powers of a/p^. The term Vsd contains spin-dependent interaction generated 
by both, perturbative and nonperturbative, potentials. Finally, the last term, Vioopin"^), comes 
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Meson 

hi 

h2 
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Mo 

AEstr 

Est 
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^SOi 
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J~'S02 

D 

1522 
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429 
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-26 

14 

15 

-13 
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4847 

675 
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5780 
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6 

7 
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Table 2: Solutions of the eigenvalue problem for the Hamiltonian (^) and the coefficients from 
equation @ for the set of parameters given in Table I. All quantities are given in MeV. 
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'Pi 

'P 2 
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1 

1 
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'Pi 


0 , 0 , - 1/4 

L 

1 

L 

0 


CO 


I/V 2 , -l/\/2, 0 

-1/2, -1/2, 1/4 


'P 2 




1/2, 1/2, 1/20 


Table 3: The matrix elements of the spin-dependent operators between P states given in the 
form: {SiL), {S 2 L), {{Sin){S 2 n)). 

from the one-loop corrections to the potential. It is given by equations (3.1) and (3.2) of the 
paper with the obvious change mi _2 hi, 2 , and we choose the renormalization scale to be 
equal to the reduced effective mass /i. Finally, to £x the Hamiltonian (^-(H), we use the values 
of the parameters listed in Table I. 

Einbein fields /ri ^2 are kept as variational parameters and the spectrum is minimized then 
with respect to them. The extremal values of /r’s play the role of the constituent masses of 
the quarks and appear dynamically due to the interaction. This feature of the given approach 
allows one to start from the current mass of the constituent (gluons also can be described in 
this formalism) and to arrive at its effective constituent mass self-consistently. However this 
simple interpretation should be considered with caution. The first source of error is neglecting 
the quark Zitterbewegung. Indeed, we neglect the negative-signed solutions for /ii ^2 expecting 
their small influence on the spectrum [Q . On the other hand, the simple quantum mechanical 
reduction of the relativistic held-theory problem given by the QCD string approach is not 
applicable for the description of chiral effects, such as the Bogoliubov-type transformation 
from bare to dressed quarks and the formation of a nontrivial chiral condensate. Therefore one 
cannot pretend to describe the pion in this framework. In realistic quantum-held-theory-based 
models each mesonic state possesses two wave functions which describe the motion forward 
and backward in time of the qq pair inside the meson Q . For the pion, which is expected to be 
strictly massless in the chiral limit, the two wave functions are of the same order of magnitude, 
so that none of them can be neglected. Luckily the backward motion is suppressed if at least 
one of the quarks is heavy p, so that one expects to arrive at reliable predictions in the case 
of heavy-light mesons. 

Since the Hamiltonian Hq, which plays the role of the zeroth order approximation for the 
problem, conserves the angular momentum L, the total spin S, and the total momentum 
J = L + S separately, then its eigenstates can be specihed as terms, with n being 

the radial quantum number. In the remainder of this letter we shall concentrate on the states 
with n = L = 1. Their masses can be represented as 

^(l2S+lp^) ^ (l2S+lp^|7f|l2S+lp_^^ ('g) 

= Mo + AEstr + EsoASiL) + ES0AS2L) + Est{ 3 {Sin){S 2 n) - {S1S2)), 
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where AEg-tr is the contribution of the string correction and the term with the spin-spin inter¬ 
action does not contribute since the wave function at the origin vanishes for orbitally excited 
states, whereas the corresponding one-loop contribution is negligible. The results of numerical 
calculations, including the values of the coefficients entering equation (H), are listed in Table 
II (see [0 for the details of the calculations). 

In Table III we give the matrix elements of the spin-tensor and spin-orbit operators between 
F-level states. Since the spin-orbit interaction mixes states with different total spin, then the 
masses of the physical states with the total momentum J = 1 are subject to a matrix equation. 




E {PPi\H\l^Pi) 
{l^Pi\H\l^Pi) -E 


= 0 . 


( 10 ) 


In Tables IV,V we give our predictions for the masses of the P-level D and B mesons and 
compare them with the predictions of other models as well as with the lattice and experimental 
data coming from various collaborations. 

From Tables IV,V one can deduce several conclusions. First, all three mentioned models 
give good description of the IP 2 states, whereas all of them fail to reproduce a very heavy 
IPo P-mesonic state reported by OPAL |^. If this experimental value is confirmed, then this 
will serve as a signal that all theoretical approaches miss something, and this question deserves 
additional careful study. On the other hand, lattice simulations give the mass 5.754GeV for this 
state [^, which is also about lOOMeV lower than the OPAL value. This stresses once again 
that the experimental situation strongly needs clarification. A similar state in the spectrum 
of D mesons is not reported yet by experimental collaborations, though all models and the 
lattice simulations give a consistent prediction for it to be around 2430 -h 2440MeV. 

Another conclusion which one can make from Tables IV,V is that there is no agreement 
concerning IPi states. Different models give different splitting patterns (see also the discussion 
in 1^). To have a better insight into the nature of this splitting let us study the heavy-quark 
limit, ttiq = mi —> 00 , analytically, which is possible in our approach. Only the coefficient 
Eso 2 , in notations of equation (pi), survives in this limit, and the expression for it reads 


a 


K 




9k^ 


lOvr/r^ LV18 


19 


+ 1e] {r ) + (r ■^ln(/ir)) 


( 11 ) 


where je = 0.5772 is the Euler constant and the averaging is performed over the zeroth-order 
wave function corresponding to both states, P 1/2 and P 3 / 2 , which are now the true 

eigenstates of the Hamiltonian [|. As discussed in |2^, || , solution of the eigenstate problem for 
the Hamiltonian (P) in this limit is given by solutions to the Schrodinger equation 


£ 


dx^ 


+ \x\ — 


A 


\x\ 


X\ = a(A),\A, 


( 12 ) 


with the reduced Coulomb-potential strength A being the solution of the equation (we put the 
light-quark current mass equal to zero for simplicity) 


A2 = (a + 2A 


da 


dX 


(13) 


^We follow the standard notations using the total momentum of the light quark, j 2 = L -I- 5'2, as the 
subscript. 
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Our 

Ref.lg 

Ref.|ll| 

Lat.|^] 

PDG 

CLEO 

IPI (.P 3 / 2 ) 

Pi 

2.428 

2.44 

2.414 

2.405 

2.422 

2.425 

IPo 

Po 

2.43 

2.4 

2.438 

2.444 



1 P 2 

P 2 

2.445 

2.5 

2.459 

2.445 

2.459 


1Pi(D/2) 

Pi 

2.469 

2.49 

2.501 

2.414 


2461 +u-U4i 

ZWi_0.034 


Table 4: Masses of the P-level D mesons in GeV. See Refs.[|!3 
values. 


1 ^ for the cited experimental 




Our 

Ref.[]^ 

Ref.|ll| 

Lat.|^] 

OPAL 

L3 

DELPHI 

CDF 

ALEPH 

1Pi(T3/2) 

Pi 

5.716 


5.719 

5.684 


5.67 


5.71 


IPo 

Po 

5.722 

5.76 

5.738 

5.754 

5.839 





1 P 2 

P 2 

5.724 

5.8 

5.733 

5.77 


5.768 

5.732 


5.739 

lPl(P/2) 

Pi 

5.741 


5.757 

5.73 

5.738 






Table 5: Masses of the P-level B mesons in GeV. See Refs. [pT 
values. 


lOf for the cited experimental 


which is Ao = 1.215 for ag = 0.39 and Aq = 1.250 for Og = 0.4. The reduced effective mass /i 
takes the value 

1 / A \ 

(r^) = (2/10-) / x^^‘^\x\{x)\ dx. (15) 

Jo 

Then the difference of the masses of the two eigenstates corresponding to j 2 = | and j 2 = | 

is 

^Pl/2 ~ ^Pz/2 = ~2^S02, (16) 

SO that the picture of the splitting depends on the sign of the coefficient (|^). Numerically this 
difference equals to +9MeV for D mesons and -|-llMeR for P’s. Of course the considered 
limit rriQ —> cx) is not realistic; it might be reasonably well justihed for the 6-quark, but not 
c-quark. In what follows the heavy-quark mass is kept hnite. 

From Tables IV,V one can see that the predictions of our method for the masses of the IPi 
states are in good agreement with the lattice calculations [Q as well as with the experimental 
data. Namely, as far as the spectrum of D mesons is concerned, we have good coincidence 
with the results of CLEO H (see Table IV). In the P-mesonic spectrum we identify the state 
Pi with the mass m(Pi) = 5.71 ± 0.02GeV, recently claimed by CDF [^, with the lightest 
member of the J = 1 doublet, whereas the heaviest one can be associated with the resonance 
reported by OPAL [|^ (see Table V). Unfortunately, experimental resolution does not enable 
one to disentangle both Pi states simultaneously, so that at present time one rather has to rely 
on available lattice simulations. In such a situation other models, as well as improved lattice 
calculations, are welcome to attack this problem to have well established and clear predictions 
for experimentalists. 

In Fig. 1 we give the splitting pattern for the IP levels for both, D and P, mesons for 
the heavy-quark mass varying from inhnity to about 1.3GeU with the vertical dashed line 
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Figure 1: Splitting pattern for the D- (left plot) and i?-mesonic (right plot) P levels as a 
function of the heavy-quark mass. The vertical dashed line corresponds to mq = rric = lAGeV 
for the left plot and tuq = ttib = 4.8Gel^ for the right one, respectively. 


giving the actual masses of c and b quarks for D and B mesons, respectively. It is also worth 
mentioning that according to our model the IPq and IP 2 levels change their ordering around 
the heavy-quark mass mg 7.9GeV for the D-like meson (left plot in Fig. 1), and around 
rriQ ~ 5.5GeV for the 5-like one (right plot in Fig. 1). 

In conclusion, let us briefly summarize the results reported in this letter. We addressed the 
question on masses and splitting pattern of the P-level D and B mesons in the method of the 
QCD string with quarks at the ends. We took into account the proper dynamics of the QCD 
string, encoded in the so-called string correction, and supplied the interquark interaction with 
the one-loop corrections adapted to the case of the self-consistently generated dynamical masses 
of the quarks. Using the standard values for the string tension, the strong coupling constant 
and the current quark masses, we calculated the spectrum of P-level D and B mesons and 
found good agreement of our results with the lattice and experimental data, including those 
reported recently. Finally, we give our predictions for the splittings between the P states. 
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